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Particle motion of a Lennard- Jones supercooled liquid near the glass transition is studied 
by molecular dynamics simulations. We analyze the wave vector dependence of relaxation 
times in the incoherent self scattering function and show that at least three different regimes 
can be identified and its scaling properties determined. The transition from one regime to 
another happens at characteristic length scales. The lengthscale associated with the onset 
of Fickian diffusion corresponds to the maximum size of heterogeneities in the system, and 
the characterisitic timescale is several times larger than the alpha relaxation time. A sec- 
ond crossover lengthscale is observed, which corresponds to the typical time and length of 
heterogeneities, in agreement with results from four point functions. The different regimes 
can be traced back in the behavior of the van Hove distribution of displacements, which 
shows a characteristic exponential regime in the heterogeneous region before the crossover 
to gaussian diffusion and should be observable in experiments. Our results show that it is 
possible to obtain characteristic length scales of heterogeneities through the computation of 
two point functions at different times. 



I. INTRODUCTION 



The basic mechanisms of diffusion in supercooled liquids and glasses are still poorly under- 
stood. Despite a considerable amount of work in the last years, different approaches to the problem 
do not always agree and raise doubts as to the ultimate relevance of the approaches themselves. 
In the landscape approacbi-^iiiSiiiSi^ the basic hypothesis is that, at low temperatures, the dynam- 
ics of a supercooled liquid or a glass is ruled by the complex topography of the potential energy 



surface, in particular by the multiplicity and complexity of local minima and saddles. The glass 
transition in this context can be viewed as the manifestation of a particular topological transition in 
the landscape, at which the mechanism of diffusion by escape through unstable directions ceases 
to be efficient and hopping through barriers becomes the dominant mechanism to get out of a re- 
stricted region in phase space. While this approach has been considerably successful and offered 
deep insights onto the nature of the dynamics in the supercooled regime, its connection with real 
space particle dynamics is not direct, and limited information on the elementary dynamical pro- 
cesses at the particle level has been obtained. At a more coarse grained level, minimal models with 
kinetic constraints have been very successful in order to get intuition and insight on elementary 
dynamical processes ^^'^^-i^'i^'i^. In this kind of models the real space dynamics can be analyzed 
in great detail and relevant mechanisms of the microscopic dynamics can be uncovered. There 
is no underlying landscape and this has raised the question of the relevance and necessity of a 
landscape based description of the dynamics of supercooled liquids. Nevertheless, in spite of their 
success in identifying many qualitative features of supercooled liquids, the extreme simplification 
in their own definition poses a limit in the description of many complex phenomena, like aging or 
thermodynamic properties of glasses. 

Mode Coupling Theor3*i^»i2iiS, is perhaps the most successful theoretical approach to the dy- 
namics of the supercooled liquid state. It predicts the existence of relevant regimes for relaxation, 
like the beta and alpha relaxations, and make quantitative predictions for temperatures above the 
glass transition. Because of the complex interplay between space and time scales in a supercooled 
liquid, most studies based on MCT have focused on time scales at fixed length scales, e.g. the anal- 
ysis of the alpha relaxation time scale is usually done focusing on the wave vector corresponding 
to the peak of the structure factor. Even the emergence of stretched exponential relaxation, a 
benchmark of glassy dynamics, is obtained from MCT only through fitting of numerical solutions 
of the complex set of equations describing time correlations. The origin of stretched exponential 
relaxation in supercooled liquids is still debated. A common scenario opposes a possibly homo- 
geneously stretched relaxation, which basically means that local regions in a supercooled liquids 
relax in the same time scale, with an heterogeneous scenario in which different regions relax with 
different characteristic times, probably exponentially, giving rise to a global stretching of the re- 
laxation. Only recently the local dynamics has begun to be accessible experimentally, through the 
emergence of new sophisticated experimental techniquesi2*22iSl and the concept of dynamical het- 
erogeneities has gained force'^^. It has been realized that complex spatio-temporal correlations 



characterize an heterogeneous system, which can be described by four point functions, correla- 
tions between two points in space at two different times22i2i2^. These dynamical susceptibilities 
show a peak at a time scale which corresponds to the typical heterogeneity in the system and this 
time scale grows when temperature is lowered towards the glass transition. An associated growing 
length-scale suggests a situation similar to a critical phenomenon. This length-scale is difficult 
to measure. Experimental results point to small or moderate values of the characteristic length of 
heterogeneities near the glass transitioni2i2^ and computer simulation results are inconclusive due 
to the limitations in the times scales reached by numerical experiments, when compared to real 
experiments^^. 

A clear signature of dynamical heterogeneity is the observation of two sets of particles with 
different degree of mobility, which allows to define "slow particles" and "fast particles" over 
particular time interval922iSii22i22i^. These two sets show up, for example, in a double peak structure 
of the distributions of particle displacements, and several indicators have been defined in order to 
locate the time at which this separation is maximal. This time scale is another way of defining the 
typical lifetime of heterogeneities. At very long time scales the supercooled liquid recovers the 
characteristics of homogeneous Brownian motion and the distribution of displacements becomes 
Gaussian. Consequently indicators of non-Gaussian behavior serve to characterize dynamical 
heterogeneity^Si^. 

Recentl}*^^, by performing a detailed study of the interplay between time and length scales, 
it was realized that it is possible to obtain characteristic length-scales of the process of diffusion 
in supercooled liquids by analyzing two point functions, namely the self part of the van Hove 
distribution function Gs{r,t) and its Fourier transform, the self incoherent scattering function 
Fs{k, t), the two most common functions used to study dynamical behavior in liquids^^. Analyzing 
the behavior of kinetically constrained lattice models, it was realized that it is possible to define 
a length-scale corresponding to the onset of Fickian diffusion in a supercooled liquid. Above this 
length the system behaves as a simple fluid and below it persistence dominates and nearly frozen 
regions are observed up to times corresponding to the alpha relaxation time, precluding the glass 
transition. This interesting observation, that stretching is dominated by persistence events up to a 
time in which Fickian diffusion sets in, can only be a rough approximation to the true behavior of 
a strongly correlated supercooled liquid, as we will see in the following. 

In this paper we show results of molecular dynamics simulations of a Lennard- Jones binary 
mixture (LJBM), focusing on the behavior of the van Hove and self scattering functions through 



a wide spectrum of time and length scale o i . Our main result is that characteristic length scales 
can be unequivocally obtained through the analysis of the wave vector dependence of relaxation 
times in two point time dependent correlation functions. We have been able to characterize at 
least three regimes: ballistic, heterogeneous and Fickian. In particular, the scale characteristic 
of the onset of Fickian diffusion is shown to correspond to the maximum size of heterogeneities. 
The other lengthscale corresponds to the typical heterogeneous behaviour and the corresponding 
time is the typical relaxation time scale. The heterogeneous regime is characterized by stretched 
behavior of the scattering function and we have observed a corresponding exponential regime in 
the van Hove distribution function of displacements. 

Our results show that typical lengthscales of heterogeneities can be obtained through analysis 
of two point functions, like the incoherent scattering function, which are easily obtained experi- 
mentally. 

In section II we introduce the model, the simulation details and show the appearance of the dif- 
ferent regimes in the inherent structures version of the van Hove distribution function. In section 
III we present our main results on the wave vector dependence of relaxation times and determina- 
tion of characteristic time and lengthscales from analysis of the self incoherent scattering function. 
In section IV we make a comparison with known experimental results. In section V we present 
some conclusions. 

n. THE LENNARD-JONES BINARY MIXTURE AND BASIC OBSERVABLES 

We performed molecular dynamics simulations on a well known Lennard- Jones binary mixture 
(LJBM) with 80 % particles of type A and 20 % particles of type B with ^aa = 1-0, ess = 0.5, 
eAB = 1-5, aAA = 1-0, (Xbb = 0.88, uab = 0.8 at a density 1.204. Most results correspond 
to systems with N=1000 particles. Some results with 130 and 10000 particles will be shown for 
discussing fine size effects. We used a cutoff radius Rc = 1.8 for the potential and obtained 
Tmct ~ 0.46 from extrapolation of diffusivity data. The simulations shown here were done at a 
single working temperature T = 0.525 = I.IATmct, at which the system was equilibrated. As 
our main aim is to analyze the diffusion dynamics of particles we map the instantaneous config- 
urations to the nearby local minima, called inherent structures^-^^. Periodically along a trajectory 
we take a configuration and let it relax with a conjugate gradient algorithm to the nearest local 
minimum. In this way we get a map of the true trajectory to a trajectory between local minima. 



As the temperature is lowered towards the glass transition, trapping in the basins of these inherent 
structures becomes important and rule the dynamics^*^. Nevertheless throughout the exposi- 
tion we make comparisons with the corresponding results from the real (instantaneous) molecular 
dynamics trajectories. 

In figure rVI we compare the behavior of the self incoherent scattering function: 



calculated from the real instantaneous configurations, with the same function calculated from the 
corresponding inherent structures. From its definition this function is wave vector dependent. In 
this figure the wave vector corresponds to the maximum of the structure factor \k\ = 7.25. The 
main difference between the two curves is the suppression of the relaxation to the plateau in the IS 
dynamics. This basically means that vibrations of the particles and rattling inside cages is strongly 
suppressed in the IS dynamics and we are left only with structural displacements. This is good for 
our purposes of looking at diffusion dynamics at low temperatures. From the instantaneous curve 
we extract as usually the a; -relaxation time scale which is approximately r ~ 120 for our system. 

The same behavior is observed in the mean squared displacement (MSD) of the particles 
R^it) = ^ Ef=i \riit) - r;(0)p shown in Fig. Note that the MSD from IS proceeds with- 
out the arrest at intermediate times. The two curves merge approximately around the a-relaxation 
time, from where the system gradually enters a normal diffusive dynamics. 

The basic quantity to characterize the displacements of particles with time is the van Hove 
distribution function^^. This function has been extensively analyzed in the context of dynamical 
heterogeneities in supercooled liquids and colloidal systems ^2iSii22i22i2fl. We consider the self part 
of the function and sum the contributions to the displacements in the three coordinate directions 
obtaining an effective one dimensional quantity defined by: 



Fsik, t) = - J2exp[ik- im - rl(0))] 



(1) 



i=l 




(2) 



(3) 



where the superscript IS means that we are considering inherent structures coordinates. The long 
time, large displacement limit corresponds to homogeneous Fickian diffusion. We were able to 
distinguish at least two regimes. In figure|V|it is shown that an exponential decay fits correctly the 



data for intermediate time scales and distances 0.5 < x < 2. We will see later when analyzing the 
wave vector dependence of relaxation times, that this regime of times and distances corresponds 
to the heterogeneous or stretched dynamical regime. An exponential decay in the distribution 
of displacements was also observed by Vogel et al.^^ and Schroeder et al.^^ when analyzing the 
displacements of particles between two consecutive inherent structures. In figure |V1 we show the 
asymptotic Gaussian regime of the van Hove function. Note that it is necessary to wait for times 
several orders longer than the alpha scale in order that almost all particles enter a Fickian regime. 
This fact was recently observed and analyzed in^. The crossover length and time scales between 
the Fickian and exponential regimes will be analyzed in grater detail in the next section. 

Recently a frozen regime was also observed in spin models of glasses with kinetic constraints^. 
In this regime a fraction of the system remains frozen up to times which increase rapidly with de- 
creasing temperature, as described by a persistence function. In the van Hove distribution the 
frozen component is reflected by a delta peak at the origin. In molecular systems a strictly frozen 
component cannot be seen due to vibrations. But from the inherent structures dynamics we ob- 
tained a delta peak for our smallest sample of N=130 particles. In this case we see that the whole 
system remains frozen for some samples up to times of the order of the typical relaxation time. 
But once a particle moves, all the others also move, although the displacements can be very small. 
As the size of the system grows, completely frozen samples became rare. For the 10000 particles 
system we have not seen the ocurrence of a frozen sample nor of a frozen region within the sam- 
ples, i.e. it is not possible to observe strictly frozen groups of particles. As the system size grows, 
the minimum possible displacement of the particles shifts continuously towards smaller distances. 
At most we must expect to see groups of particles with distinctive mobility, i.e. slower and faster 
particles. 

m. THE SELF SCATTERING FUNCTION AND CHARACTERISTIC TIME-LENGTH 

SCALES OF HETEROGENEITIES 

Assuming isotropy of space we computed the following one dimensional self incoherent scat- 
tering function: 

1 ^ 

FsHk,t) = -Eexp[.A;(xf (t) -xf (0))] (4) 
i=i 

in which the superscript IS means that the function is calculated from inherent structures configu- 
rations. 



In figure |V1 we show four selected curves corresponding to different wave vectors in order to 
illustrate the different regimes of relaxation. For small k the relaxation is exponential correspond- 
ing to Fickian diffusion of these modes. For wave vectors around the peak of the structure factor 
(/c = 7 is shown in the figure) or larger, the long time decay can be well fitted by a stretched 
exponential with a k dependent exponent which will be analyzed below. The regime of large A;'s 
correponds to the system being around the basin of a single inherent structure, and the decay can 
be fitted with two stretched exponentials, one for the short time and another for the long time 
regimes of the relaxation. This analysis was performed and discussed in detail in^Siit. In figure IVI 
we show the main result of this work, the dependence of the relaxation time r(A;) on wave vector. 
Relaxation times for each k were obtained as usual, as the time at which the self scattering func- 
tion decays to 1/e of its initial value. In this plot we compare the results from inherent structures 
(IS) and real coordinates for three system sizes: = 130, = 1000 and = 10000. We see 
that the data for the two larger systems almost coincide, so one can be confident that already for 

= 1000 there are no finite size effects in the scaling regimes. This figure is very interesting^, 
showing different scaling regimes between time and length scales. For each system size and for the 
smallest wave vectors, the results from IS and real coordinates coincide, while they differ strongly 
at large wave vectors, as expected. 

For wave vectors k < 2 and times r > 1000, a Gaussian scaling r ^ 1/A;^ is observed: this is 
the Fickian diffusion regime (see the inset in figure IVJ. At A; ^ 2 a crossover to another scaling 
form is observed: the system enters an anomalous diffusion regime with scaling r ~ 1/ k^-^ . This 
regime extends between 2 < A; < 10 and times between 50 < r < 1000. This is the regime 
in which the dynamics is heterogeneous and the relaxation of correlation functions is stretched. 
Comparing the length and time windows in this regime we realize that they correspond to the 
exponential decay of the van Hove distribution (see figure IVT). 

For A; > 10 the relaxation times decay rapidly in a narrow interval of wave vectors. Note that 
this effect of rapid decay is much more pronounced in the data corresponding to the real coordi- 
nates than in the inherent structures data. This is due to the cage effect, in which the diffusion 
is temporarily halted, seen in the figure as a narrow interval in k values. For scales larger than 
A; > 20 (data not shown), a linear scaling works well: the diffusion is ballistic. While the ballistic 
regime is nicely observed in the real dynamics, the corresponding curves for the inherent struc- 
tures dynamics show instead a slow decay of relaxation times with k, because the fast motion of 
the particles is filtered out. The second crossover, between the subdiffusion and the faster decay 



happens at a time of the order of the alpha relaxation time. This time corresponds to the peak in the 
four point susceptibility and the corresponding length scale to the typical size of heterogeneities. 

The crossover between Fickian and anomalous diffusion at nearly kp — 2 allows to define a 
characteristic length scale /i? ~ 3 as the scale for the onset of Fickian diffusion at the temperature 
considered. Below this scale cooperativity rules the motion of individual particles and the dynam- 
ics is heterogeneous. This length scale corresponds to the maximum size of heterogeneities in the 
system. The departure from anomalous diffusion at nearly — 10 allows to define a second 
length scale ~ 0.6, as the typical size of heterogenous regions. Note that this size is very 
small at this temperature. It would be extremely interesting to repeat this analysis for different 
temperatures approaching and crossing the mode coupling transition temperature. 

In figure IVl we show the dependence of the stretching exponent P^^{k) on wave vector. It 
decreases monotonously from exponential (/5 = 1, Fickian) behavior at small k to an apparent 
saturation around /? ~ 0.3 at large k. Due to the complex t and k dependence of the stretching 
exponent, we were not able to collapse the data of the self scattering function for different wave 
vectors onto a single master curve. 

rV. EXPERIMENTAL EVIDENCE 

There is a relatively large literature reporting results on heterogeneous dynamics in glasses and 
the search for characteristic time and length scales^^. Nevertheless there are, to our knowledge, 
only a few papers reporting detailed measurements of particle displacements in supercooled liquids 
and glasses^SiSi'^. 

In figure 13 we show that in the real dynamics, for a time less than the alpha relaxation time, 
the van Hove function can be fitted by a gaussian contribution from vibrations at small distances 
plus an exponential contribution at large distances. This figure can be compared, e.g. with figure 
3 of^\ where experimental data from colloidal supercooled fluids and glasses was fitted with an 
stretched exponential with exponent ~ 0.8, although an exponential fit can probably work also in 
that case. 

In another, similar experiment, Kegel et aP°, using time-resolved fluorescence confocal scan- 
ning microscopy, also measured the self part of the van Hove function. They classified the particles 
in two subsets: one fast and one slow, and fitted the data for both subsets with two Gaussians. The 
fits are very good, as recently observed for the same time scale in simulations of the LJBM^. The 



difference in the behavior of the distribution function between both experimental results probably 
reflects the different time regimes in which the measurements were done in each case. It would 
be very interesting to have accurate measurements in the whole time span between the beta re- 
laxation scales up to several times the alpha scale in order to get from that an estimation of the 
characteristic length scales discussed in this work. 



V. CONCLUSION AND PERSPECTIVES 



The mechanisms of particle diffusion in deeply supercooled liquids are still poorly understood. 
In particular, the region of intermediate time and length scales, in which diffusion is anomalous, is 
still waiting for a complete theoretical description. In the meantime, new highly precise measure- 
ments probing local dynamics and extensive computer simulations are giving important insights 
on the basic mechanisms that underlie particle motion. By means of molecular dynamics simu- 
lation on a Lennard- Jones supercooled liquid we showed that the van Hove distribution and self 
scattering functions still bring us new and rich information. They show at least three well de- 
fined regimes on different time and length scales. On very short times and lengths the motion 
of particles is ballistic. Then heterogeneities develop and at times of the order of the relaxation 
time a well defined crossover length In can be obtained from the k dependence of the relaxation 
times. This length corresponds to the typical size of spatial heterogeneities. In this region the 
diffusion is anomalous and time correlations decay in a stretched exponential way. Also in this 
space-time regime the van Hove distribution shows a well defined exponential decay. Up to our 
knowledge this exponential decay has not been observed nor analyzed in other models, like the 
much studied kinetically constrained lattice models, but is clearly present in experimental results 
on colloidal systems. At very long time scales, several orders larger than the alpha relaxation 
scale, Fickian diffusion sets in, and the distribution of displacements slowly converges to a Gaus- 
sian. The crossover from heterogeneous to Fickian dynamics is also clearly observed in the k 
dependence of the self scattering function, and a second typical length Ip can be obtained. Scaling 
forms typical of normal and anomalous diffusion can be seen for distances larger and smaller than 
If, respectively. An intersting study on the temperature dependence of the onset time for Fickian 
diffusion was recently done by Szamel et al.-^. In that work the onset time for Fickian diffusion is 
defined as the time at which a fixed deviation from gaussianity was observed in the distribution of 
the logarithm of particle displacements P(logio5r,t). Its temperature dependence was analyzed 



and compared with the corresponding behavior of other characteristic time scales. A particularly 
interesting result is the observation that the ratio between the Fickian diffusion onset time Tp and 
the alpha relaxation time growths with decreasing temperature, but tends to saturate near the 
Mode Coupling transition temperature Tc. This may indicate that the mechanisms behind both 
time scales are essentially the same as the crossover temperature is approached. 

Although the identification of characteristic dynamical lengths can be naturally introduced 
through four point correlation and response functions, we have shown that the relevant infor- 
mation can already be obtained from two point functions, like the van Hove distribution and self 
scattering correlation functions. To understand the emergence of the different regimes observed 
in the simulations from a single microscopic model, or unified theoretical framework, is still a big 
challenge. 
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Captions for figures 

• Figure 1: Self incoherent scattering function for the LJBM at T=0.525 for a wave vector 
corresponding to the peak of the structure factor k=7.25, from instantaneous and inherent 
structures coordinates. 

• Figure 2: Mean squared displacement at T=0.525 for instantaneous and inherent structures 
coordinates. 

• Figure 3: The van Hove distribution in the exponential regime. Continuous lines are expo- 
nential fits. 

• Figure 4: The van Hove distribution in the Fickian regime. Continuous lines are gaussian 
fits. 

• Figure 5: The self incoherent scattering function for four different wave vectors characteris- 
tic of different scaling regimes. In the legend the corresponding fitting functions (continuous 
lines). 

• Figure 6: Wave vector dependence of relaxation times for three different sizes, inherent and 
instantaneous dynamics. Different scaling regimes indicated with full lines. Inset: zoom of 
the Fickian crossover region. 

• Figure 7: Wave vector dependence of the stretching exponent p. 

• Figure 8: The van Hove distribution from instantaneous coordinates att = 32. The data can 
be fitted by a Gaussian at small distances plus an exponential decay at larger distances. 
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